Weakly coupled s = 1/2 quantum spin singlets in Ba 3 Cr 2 
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Using single crystal inelastic neutron scattering with and without application of an external mag- 
netic field and powder neutron diffraction, we have characterized magnetic interactions in BaaCi^Os- 
Even without field, we found that there exist three singlet-to-triplet excitation modes in (h, h, I) scat- 
tering plane. Our complete analysis shows that the three modes are due to spatially anisotropic 
interdimer interactions that are induced by local distortions of the tetrahedron of oxygens sur- 
rounding the Jahn-Teller active Cr 5+ (3d 1 ). The strong intradimer coupling of Jo = 2.38(2) meV 
and weak interdimer interactions (|Ji n t cr | < 0.52(2) meV) makes BaaC^Os a good model system 
for weakly-coupled s = 1/2 quantum spin dimers. 



PACS numbers: 75.10.Jm 

For the past two decades, the field of quantum mag- 
nets has provided a number of new exotic collective phe- 
nomena [I]. Progress has always been moderated by the 
availability of suitable spin 1/2 materials with distinct 
structural motifs. The ubiquitous Cu 2+ (3d 9 , s = 1/2) ion 
is often relied upon as the platform for new materials de- 
sign. Magnetic field-induced quantum phase transitions 
are a typical example. The field-induced condensation 
of magnons have been experimentally observed in cou- 
pled quantum (s = 1/2) dimer systems based on Cu 2+ 
ions, such as TlCuCl 3 HUSH] and BaCuSi 2 6 [5], which 
are adequately described by the Bose-Einstein conden- 
sation (BEC) theory [21 [HIE]. However, the robustness 
of such descriptions can only be truly evaluated with in- 
vestigation into complementary materials. This calls for 
finding other good model systems of quantum (s = 1/2) 
dimers, especially based on non-Cu 2+ ions. 

Recently, a new class of spin dimer systems has been 
reported with the general formula, A 3 B 2 8 [51 151 ITD1 ITT] . 
where A is an alkaline earth metal such as Ba 2+ or Sr 2+ 
and B is a 5+ transition metal ion that forms dimers 
along the c-axis (see Fig. [l|b)). Interestingly, a frus- 
trating triangular lattice is formed by the dimers in the 
a&-plane (see Fig.HTa)) [5]. For Ba3Mn 2 0s with orbitally 
nondegenerate Mn° + (3d 2 ,s = 1), specific heat measure- 
ments under H provide experimental evidence of a con- 
densation of magnons [D] ■ In the field-induced antiferro- 
magnetic (AFM) state, a magnetization plateau phase 
was observed as a function of H [8J, typical for frus- 
trated magnets [12] . A recent inelastic neutron scatter- 
ing study showed that the magnetic interactions between 
dimers in Ba3Mn 2 0s are effectively two-dimensional in 
nature as a result of this interplane frustration |13j . More 
recently, the related chromate (BajSr^Q^Og with the 
quantum (s = 1/2) spin from the orbitally degenerate 
CrS+^d^s = 1/2) ion have been made [10l [11]. These 
compositions do not exhibit long range order above 1.5 K 
and their bulk susceptibility data at low temperatures 



observes the quantum dimer model. The application of 
an external magnetic field at 1.6 K rapidly increases the 
bulk magnetization in Ba 3 Cr 2 08 above H cl ~ 12 T to 
its saturation value at H c \ « 23 T [ID]. In contrast to 
the manganate, no magnetization plateau was observed. 

Here we report our inelastic neutron scattering data 
obtained from single crystals of Ba3Cr 2 0s with and with- 
out application of an external magnetic field. We have 
found that even without field there are three singlet-to- 
triplet excitation modes in (h, h, I) directions, which con- 
trasts with the single mode observed in Ba3Mn 2 0s- Our 




FIG. 1: Schematic diagram of the crystal structure of 
Ba3Cr2 08 that shows (a) the triangular network of Cr + 
dimers in the afc-plane and their relative stacking along the 
c-axis in an ABCABC arrangement, (b) Two oriented Cr 
tetrahedra bound with four oxygen (O 2- ) ions to afford the 
Cr 5+ dimer arrangement with magnetic exchange, Jo, which 
is at a distance of do = 3.97A. (c) The ab-plane projection of 
the Cr 5+ ions, showing the alternative exchange pathways. 



detailed analysis shows that the three modes come from 
spatially anisotropic interdimer interactions in three dif- 
ferent crystal domains. We have also performed neutron 
diffraction measurements on a powder sample and found 
experimental evidence of local lattice distortion that is 
due to the orbital degeneracy of the Cr 5+ (3d 1 ) ion. We 
consider that the local distortion is the origin of the spa- 
tially anisotropic interdimer interactions. Our complete 
analysis of dispersions of the singlet-to-triplet excitations 
show that this system is an excellent model system for 
three-dimensional weakly coupled dimers with dominant 
intradimer coupling (Jo — 2.38(2) meV) and weak inter- 
dimer couplings less than 0.52(2) meV. Thus, this system 
provides a new testing ground for theories of the conden- 
sation of magnons. 

Single crystals of Ba 3 Cr 2 0s were grown by the floating 
zone method using a four mirror image furnace equipped 
with halogen lapms. Since Ba3Cr 2 0s is congruently 
melting, sintered stoichiometric rods were used. We have 
performed inelastic neutron scattering measurements on 
single crystals with a total weight of 450 mg at the cold- 
neutron triple-axis spectrometer (SPINS) and neutron 
diffraction measurements on a 15 g powder sample of 
Ba3Cr 2 0g at the BT1 neutron powder diffractometer, 
located at the National Institute of Standards and Tech- 
nology (NIST) Center for Neutron Research. At SPINS, 
the crystals were mounted in the (hhl) scattering plane. 
A horizontally focusing analyzer mode was set at the fi- 
nal energy of Ej = 5 meV and to cover the scattering 
angle of about 5 degrees at a time. This increased the 
sensitivity of the instrument and allowed us to map out 
magnon dispersions along high symmetry directions. A 
liquid nitrogen cooled Be filter was placed before the an- 
alyzer to eliminate high order neutron contaminations. 
At SPINS, the samples were cooled in a liquid 4 He cryo- 
stat for the measurement without field, and in a 11.5 T 
superconducting magnet for the field measurements. At 
BT1, a 4 K closed-cycle refrigerator (displex) was used. 

In order to investigate the intra- and inter-dimer in- 
teractions in Ba3Cr 2 0s, we have performed inelastic 
neutron scattering measurements on single crystals at 
1.7 K. Several constant Q-scans were performed to map 
the dispersion of the singlet-to-triplet excitations. Fig- 
ure [2] shows the results along three high symmetry direc- 
tions. Along the (0,0, 1) direction, a nearly flat excita- 
tion mode centered ~ 2 meV is observed with a band- 
wdith of ~ 0.7 meV (see Fig. [2ja) ) . This dispersion is 
due to the weak interdimer interactions, J\ and J4, be- 
tween the dimer planes (see Fig. [lj. The (ft, ft, 3) and 
(ft, ft, 3ft) directions, that probe the inplane interdimer 
interactions, J 2 , J3, as well as the out-of-plane ones 
show three unexpected excitation modes (see Fig. [2|b) 
and (c)). In contrast, Ba3Mn 2 0s displays only a sin- 
gle triplet excitations along all the symmetric directions 
measured [13]. Magnetic excitations from weakly cou- 
pled dimer systems can be well explained within the 
random phase approximation (RPA) that gives the dis- 
persion relation [H [H [TS], ftw(Q) ~ y/J$ + J l(Q), 
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FIG. 2: Contour maps of dispersion relations along 
(a) (0,0,/), (b) (ft, ft, 3) and (c) (ft, ft, 3ft) directions. The 
data were taken by performing constant-Q scans of BaaC^Os 
single crystals at 1.7 K. The point symbols represent the 
peak positions obtained by fitting the data with simple Gaus- 
sians. (d)-(e) Calculated dispersions and intensities based on 
a model of spatially anisotropic interdimer interactions, which 
is described in the text. 



where j(Q) = J2i J{Ri)e~ lC} ' Ri and J(Ri) represents in- 
terdimer interactions. If the interdimer interactions are 
spatially isotropic, as shown in Fig. [lja), then the four 
nearest interactions, J%, J 2 , J3, J4, gives 
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such that without field only a single triplet excitation 
along any of (ft, ft, I) directions is expected, as recently 
observed in the related s = 1 system Ba3Mn 2 0s [T3"] . 

This poses the question as to the origin of the complex 
excitations observed in this present study on Ba 3 Cr 2 08. 
One possibility would be that the triply-degenerate 
singlet-to-triplet excitation mode is already split into 
three modes in the (ft, ft, I) scans as a result of single 
anisotropy terms. If this is the case then each excita- 
tion in the (ft, ft, I) scan are not degenerate and therefore 
remain single under an application of an external mag- 
netic field. This conjecture was tested by performing fur- 
ther inelastic neutron experiments at constant-Q under 
an external magnetic field, which are shown in Fig. [3] In 
Fig. ga), the single mode at Q = (0, 0, 4.5) at H = T 
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FIG. 3: Constant Q scan with several different magnetic 
fields (H) (a) at Q = (0,0,4.5) and (b) at Q = (0.5,0.5,3) 
scan, and (c) ff-dependence of the peak positions in tko (c) 
at Q = (0,0,4.5) and (d) at Q = (0.5,0.5,3). The lines 
represent the Zeeman splitting. 



splits into three modes for finite values of field. As H 
increases, the Zeeman energy splitting increases as ex- 
pected. Furthermore, Fig. [3jb) shows the two strong 
modes at Q = (0.5, 0.5, 3) that appeared at Hlo =1.4 meV 
and 3.05 meV without field also revert to three modes 
when H is applied. This demonstrates unambiguously 
that both modes are triply degenerate singlet-to-triplet 
excitations. 

Another possibility is that in analogy to the acoustic 
and optical vibrational modes that exist in a lattice with 
more than one ion in a unit cell, the origin might be 
an acoustic and two optical magnon modes as a result 
of the three bi-layers in the chemical unit cell. Such an 
optical magnon mode has been observed in Cs 3 Cr 2 Brg, 
where Cr 3+ (3d 3 ;s = 3/2) ions form two bi-layer trian- 
gular planes [H]. The optical mode, however, appeared 
along certain symmetry directions with much weaker in- 
tensity than the acoustic mode. For Cs 3 Cr 2 Brg, Leuen- 
berger et al. have calculated, using a Green's function 
technique, the dispersions and intensities of the acoustic 
and optical modes. They showed that the intensity of the 
acoustic mode is proportional to 1 + cos(p ■ t + <f>) while 
that of the optical mode is related to 1 — cos(p ■ r + <j>) 
where p is the displacement vector between two different 
sublattices (planes) , t a reciprocal lattice vector, and <j> is 
the phase of the Fourier sum of the interplane interactions 
which is complex in the case of Cs 3 Cr 2 Brg [TJ]. Their cal- 
culations were consistent with their data. We have im- 
plemented their model for the present case, Ba 3 Cr 2 8 , 
which suggests the presence of one acoustic and two opti- 
cal modes. However, for this case <f> turns out to be zero 
because of the inversion symmetry at the center of the 
dimer and therefore the Fourier sum of the interplane in- 
teractions becomes real. Since r(h,k,l) must satisfy the 



reflection condition — h + k + I = 'in with an integer n 
for the B2m symmetry and p = (§, |, |), cos(p ■ r) is 
always 1. Therefore, the optical modes in Ba 3 Cr 2 08 are 
not present, and the observed excitations are all acous- 
tic in nature, which is consistent with the observation 
of only one single triplet mode in the related compound 
Ba 3 Mn 2 8 [13]. 

A third possibility is that the interdimer interactions 
are not spatially isotropic and the three modes arise from 
three different crystal domains. To test this scenario, 
we have considered a set of spatially anisotropic inter- 
dimer interactions, J 1 ,J 1 ,J 1 , J 2 , J 2 , J 2 > ^4> J 4 > J 4 > as 
depicted in Fig. [TJc). J 3 is neglected because we cannot 
observe J 3 separately but we only observe J 2 — J 3 . Then, 
for one crystal domain 
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produce the same 
dispersion along the (0,0,1) direction, whereas along the 
(h, h, 3) and (h, h, 3h) directions each domain will yield a 
unique dispersion as shown in different lines in Fig. |2^b) 
and (c). Fig. [2^d)-(f) are calculated dispersion relations 
and demonstrate that the experimental data can be sim- 
ulated by this model when J = 2.38(1), j[ = 0.08(2), 
J'[ = -0.15(2), J" = 0.10(1), J' 2 = 0.10(2), 
J 2 = 0.07(2), J 2 " = -0.52(2), J 4 = 0.04(2), 
j'l = 0.10(2), J 4 = 0.09(2). All values are given in meV. 

A central question arises as to what is the origin of the 
spatially anisotropic interdimer exchange interactions in 
Ba 3 Cr 2 C>8? In order to address this issue, we have per- 
formed neutron diffraction measurements on a powder 
sample at various temperatures from 4 K to 300 K. We 
found that there is no evidence of a reduction of the am- 
bient i?3m crystal structure symmetry on cooling down 
to 4 K. The displacement parameters, By, of the apical 
oxygen, however, show anomalous behavior as a func- 
tion of T. Firstly, a significant improvement in the Ri- 
etveld fit to the diffraction pattern is obtained with the 
use of anisotropic values for the displacement parameters, 
Bij, the temperature dependence of which are shown in 
Fig. |4];b). Secondly, upon cooling the ab-components of 
the displacement parameter, B±i = B22, for the apical 
(6c) oxygen decreases down to ~ 100 K but then increases 
on further cooling (see red circles in Fig.[4jb)) while the 
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FIG. 4: (a) Observed, calculated and difference intensities as 
a function of two-theta as obtained from Rietveld refinement 
of the powder neutron diffraction data at 4 K and 300 K. 
Expected reflections are represented by green tickmarks. (b) 
T-dependence of the anisotropic displacement parameter, Bij, 
for the apical (red circles) and for the planar (black circles) 
oxygens. 



c-component , B33, decreases below 100 K. In contrast, 
for the planar (18h) oxygen shows a very weak In- 
dependence below 100 K. This suggests that the apical 
oxygen that is located at a high symmetric 6c(0, 0, z) 
site locally moves in the ab-plane for T < 100 K. The 
local disorder is presumably a result of the e g orbital 
degeneracy within the Cr 5+ (3rf 1 ) tetrahedral ion [15] , 
which has a tendency to affect the local oxygen envi- 
ronment and thereby removing the degeneracy. Even 
though the Jahn- Teller distortion does not seem to lower 
the bulk crystal symmetry in BasC^Os, the anomaly 
in Bij indicates the presence of local disorder, which in- 
duces the spatial exchange anisotropy that was observed 
in our inelastic neutron scattering measurements. Very 
recently, Chapon et al. |17j have studied a related ma- 
terial Sr3Cr 2 0s using a powder neutron diffraction tech- 
nique and found that at 1.6 K the crystal structure is 
distorted to a monoclinic C2/c symmetry that also in- 
duces similar spatial exchange anisotropy, which supports 
our scenario of the local disorder in BaaC^OV We ex- 
pect that S^C^Os would also exhibit more than one 
single triplet excitation mode without field, which has 
very recently been confirmed by inelastic neutron scat- 
tering measurements on a single crystal of Sr^CroO«[l8]. 
In comparison, BaaMnaOg lacks any orbital degeneracy, 



and thus has spatially isotropic Js that leads to the single 
triplet excitations as observed |13j . 

In summary, the triplet excitations in the spin gap sys- 
tem, Ba3Cr20§, are measured by inelastic neutron scat- 
tering on single crystals while its crystal structure is de- 
termined from powder neutron diffraction experiments. 
The combination of these results show that Ba3Cr 2 0g 

TABLE I: The refined crystal structural parameters of 
BasCr208 derived from the data shown in Fig. |4ja) using the 
program, Fullprof. The symmetry was found to be R3m at 
all temperatures. W represents the crystallographic Wyck- 
off position, \ 2 = {/°ks,i _ Ical.i} 2 /A/o bSjl and R F = 

100 x J2 IW»)M) - / ca i(M,0l/]T \I ohs (h,k,l)\ where 

h,k,l h,k,l 

7obs and 7 ca i are the observed and calculated intensity, respec- 
tively. 



Atom(TT) x 



B ia o(A 2 ) 



300 K, x = 1-30, Rp2 
a = b = 5.74171(6)A, c 




0.82832(10) 




Bal(3a) 
Ba2(6c) 
Cr(6c) 
01(18/i) 
02(6c) 



= 2.99 

= 21.38736(30)1 



0.20587(8) 

0.40704(8) 

0.17168(10) 0.89874(3) 

0.32847(6) 



1.015(36) 
0.698(26) 
0.689(37) 



4 K, x = 2.86, R F 2 = 
a = b = 5.72077(8)1, c 



Bal(3a) 

Ba2(6c) 

Cr(6c) 

01(18/i) 

02(6c) 







0.82831(13) 




4.49 

= 21.35857(37)1 



0.20576(11) 

0.40648(11) 

0.17169(13) 0.89866(4) 

0.32842(8) 



0.466(45) 
0.216(33) 
0.110(47) 



Anisotropic displacement parameter (Bij x 10 s A 2 ) 


Atom(W) 


■Bll B22 Byi 


B\2 B\z B23 


01(18/i) 
02(6c) 


300 K 

1242(28) 1242(28) 62(2) 
2113(32) 2113(32) 22(3) 


945(31) 12(1) -12(1) 
-1056(16) 


01(18/0 
02 (6c) 


4 K 

475(36) 475(36) 31(2) 
1767(42) 1767(42) 12(4) 


495(41) 1(5) -1(5) 
-884(21) 



is a good model system for weakly-coupled s = 1/2 
quantum spin dimers with the intradimer coupling of 
JO = 2.38(2) meV and the interdimer couplings less than 
0.52(2) meV, which are spatially anisotropic due to local 
disorder induced by the Jahn- Teller active Cr 5+ ions. 
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